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In a recent papei*^ we have shown that the dissipation driven quantum phase transition of the 
2D xy model represents a universality class where the correlations at criticality is local in space 
and power law in time. Here we provide a detailed analysis of the model. The local criticality is 
brought about by the decoupling of infrared singularities in space and time. The former leads to a 
Kosterlitz Thouless transition whereby the excitations of the transverse component of the velocity 
field (vortices) unbind in space. The latter on the other hand leads to a transition among excitations 
(warps) in the longitudinal component of the velocity field, which unbind in time. The quantum 
Ashkin- Teller model, with which the observed loop order in the Cuprates is described maps in the 
critical regime to the quantum xy model. We also discuss other models which are expected to have 
similar properties. 



The dissipative quantum 2D xy model was 
introduced 2,3 to describe experiments on ultrathin 
granular Superconducting films, where it was observed 
that above a normal-state sheet resistance of order 
Rq = h/Ae 2 the resistivity does not decrease towards 
even at the lowest temperatures studied^.. The 
granular superconductor is represented by an array of 
islands with finite superfluid density with Josephson 
coupling among the nearest neighbors. Dissipation of the 
Caldeira-Leggetl£ form is included, whereby the phase 
difference between the nearest neighbors is coupled to an 
external bosonic bath, the external bath effectively act 
as a resistor, R, shunting nearest neighbor islands and 
energy is dissipated at the rate V 2 / R, where V is the 
voltage drop which is proportional to the rate of change 
of the phase difference between the grains. A quantum 
phase transition occurs as the value of the resistance is 
tuned. For small values of the resistance, phase slips, 
whereby a voltage develops across the junction between 
two grains, cost too much energy and so lead to an 
ordered phase. On the other hand, for large value of the 
resistance, phase slip events are energetically favored 
and global superconductivity is destroyed. 

In the absence of dissipation, the quantum model at 
T — > maps to a 3d xy model. It was suggested that 
the physics in the presence of dissipation may not be 
continuously connected to the 3d xy universality class^. 
It was also suggested that there is a possibility that the 
singularity introduced by dissipation may be linked to 
a transition driven by the unbinding in time of phase 
slips 7 . The calculations are done in the weak coupling 
regime where one considers the Josephson coupling per- 
turbatively. The same approach was previously used to 
study the phase transition and correlation functions 2 -^. 
Arguments were given also that the critical fluctuations 
due to the unbinding of the phase slips in time which 
have long-time singularities may be local in spaced. How- 
ever later work by some of the same authors presents a 
rather different pictured. Moreover none of these investi- 
gations into the dissipative transition are applicable near 
the critical point. 



In a previous paper we have shown that the the quan- 
tum phase transition is obtained due to the prolifera- 
tion of a new kind of excitation, which we termed warps, 
whose physics is indeed related to the phase slips in time; 
more specifically they represent events in time where 
a change in the longitudinal component of the velocity 
fields occur. The possibility that the longitudinal com- 
ponent may lead to the introduction of sources and sinks 
which control the critical behavior was also recently em- 
phasized in the context of one dimensional system o 10 ' 11 . 
The correct identification of the relevant degrees of free- 
dom in higher dimensions turns out to be interesting. 
The introduction of warps simplifies the analysis of the 
problem so that the correlation function of the order pa- 
rameter can be obtained near the critical point. We find 
spatially local correlations with an infrared divergence 
(oc 1/t) in temporal correlations at the critical point. 
We also the cross-over from the quantum critical to the 
quantum- regime. 

The motivation for this study is that the existence of 
such a critical point naturally explains some of the most 
important features of the phase diagram of the Cuprate 
superconductors. Polarized neutron scatterin g 12 i 13 ' 14 i 15 
in three distinct family of Cuprates have confirmed previ- 
ous finding in a third family of Cuprates through dichro- 
ism in Angle Resolved Photoemission spectral, which 
reported the observation of time reversal violation in the 
pseudogapped phase of underdoped cuprates. In par- 
ticular, the experiments lend support to the idea that 
the ground state in the pseudogapped regime is one 
with ordered current loops without loss of translation 
symmetr y 1 7 1 1 8 ' 1 9 . The low energy sector of this theory 
contains four states that are related by time reversal and 
reflection along the x or y axis. Alternatively, the four 
states can be represented by two Ising variable in each 
unit cell, the ± values representing currents in the x and 
y directions. A classical statistical mechanical model for 
such degrees of freedom that supports a phase transition 
without any divergence in the specific hea t 20 ' 21 is the 
Ashkin- Teller model 2 ^. Such a model is known to un- 
dergo a Gaussian phase transition 2 ^. A quantum gener- 
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alization of the classical model in terms of operators that 
induce fluctuations within the low energy phase space of 
the four states, leads to the dissipative 2D xy model stud- 
ied here. The local criticality with a 1/t singularity is the 
fluctuation spectrum necessary to explain the anomalous 
properties in a number of properties such a resistivity, 
nuclear relaxation rate and optical conductivity in the 
normal state near optimal doping 2 ^ abutting the pseu- 
dogapped phase. At long wavelengths, such correlations 
are directly observed by Raman scattering experiments. 
Recently, we have also shown that such fluctuations cou- 
ple to fermions to promote pairing in the <i-symmetry 
channel. 25 

It is worth pointing out that a number of classical sta- 
tistical mechanical models such as the 6-vertex and 8- 
vertex models also have Gaussian criticality in part of 
their phase diagram, as well as phase transitions which 
do not manifest any anomaly in their specific heat. Such 
models when augmented with ohmic dissipation do be- 
long to the same universality class considered here. Lo- 
cal quantum criticality happens also to be a hallmark 
of the quantum phase transitions in heavy fermion o 26 ' 27 , 
although it is not clear to us how the class of models 
solved here may be connected to the microscopic models 
for heavy fermions. 

This paper is organized as follows: The dissipative 
quantum XY model and the Villain transformation are 
discussed in sec II. In section III, we show how dissi- 
pation changes the critical properties of the model. We 
show that the action can be written in a quadratic form in 
terms of two orthogonal topological excitations, vortices 
which interact only in space and warps which interact 
only in time. We also discuss the difficulties in analyzing 
the model based on other choice of variables. In sec IV 
we discuss the properties of the warps. Sec V deals with 
the scaling equations for the vortices and warps and the 
phase transition which occurs through the proliferation 
of warps. In sec VI we present a detailed calculation of 
the fluctuation spectra near the quantum critical point. 



The spectra has the same form as that proposed to ex- 
plain the normal state marginal fermi-liquid properties 
of the high T c cuprates. In sec VII we present the zero 
temperature phase diagram. The connection to Ashkin 
Teller model and cuprates is presented in sec VIII. The 
effect of fourfold anisotropy on the critical properties is 
anayzed in sec IX. In sec X we show how dissipation 
arises in cuprates. In sec XI we discuss other forms of 
dissipation and how they may lead to different forms of 
singularities in the fluctuation spectra. 

I. DISSIPATIVE QUANTUM 2DXY MODEL 

The classical 2d xy model consists of U(l) degrees of 
freedom, represented by an angle 8, living on the sites 
of a regular lattice, assumed to be a square lattice here, 
with a nearest neighbor interaction of the form 

H= (1) 

m 

where J is Josephson coupling. Since a continu- 
ous symmetry cannot be spontaneously broken in two 
dimension s 28 ' 29 , this model does not support a long range 
ordered phase. Nevertheless a phase transition does oc- 
cur at finite temperature where the correlation function 
of the order parameter e 10 changes from exponential to 
power law. This is the Kosterlitz-Thouless-Berezinskii 
transitio n 30 ' 31 . The physics of this phase transition is 
better understood in terms of the topological defects of 
the system. To do so we follow the standard procedure of 
using the Villain transform and integrating out the phase 
degrees of freedom. We include the algebra here so as to 
make the later discussion of the quantum version easier 
to follow. 

The Villain transformation expresses the periodic func- 
tion in terms of a periodic Gaussian: 
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where rriij are integers that live on the links of the original 
square lattice. On a square lattice we choose to identify 
each site as i = x, y where x and y are integers. We can 
combine the two link variables m x . y:x +i^ y and m x . y;xy +i 
into one two component vector m 1)S that lives on the site 
x, y of the lattice. The choice is purely for convenience 
and none of the results depend on how one chooses to 
organize the degrees of freedom of the model. Instead of 
introducing a new field to linearize the quadratic term via 



the Hubbard-Stratonovich transformation, we follow an 
alternative procedure. We expand the quadratic term 
and transform to Fourier space. Keeping the leading 
quadratic term 8 XtV — Q x +\^ y ~ — a\7 x 6 xy , where a is the 
lattice constant we get 



(0 xy ~ Ox+i,y ~ 27rm* J w a 2 X7 2 J xy + AnaV x 6 xy m x xy 

+ in 2 m^ y , (3) 
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where m x is the x component of the vector field given 
by the integer y. x +i y . Fourier transforming and in- 
tegrating out the field 9^ in terms of which the action is 
quadratic we get 
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Combining the two terms and Fourier transforming to 
real space we get 



In this model, a phase transition occurs from a state at 
high temperatures where the vortices are free to a low 
temperature state where they are bound in pairs. 

We now turn to the description of the Quantum ver- 
sion of the 2d xy model. First we consider the model in 
the absence of dissipation. The quantum 2d xy model 
is written in terms of operators 9 and its conjugate, the 
number operator n. The Hamiltonian is of the form 
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Since the vector field is two dimensional the 
scalar and is the vorticity of the vector field, 
integer charges V x m r = p r we get the standard action 
of the two dimensional coloumb gas. On the lattice this 
definition is equivalent to 
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= iSij. Since 



where C is the capacitance and [nt, 
rii = id/d8i, the first term is the kinetic energy of the 
phase with a mass C . There is competition between the 
kinetic energy and the potential energy terms, the former 
minimized by a state where 9i is disordered stabilizing 
an insulating phase while the latter minimized by a fixed 
value of 9i stabilizing a superconductor. The partition 
function can be recast in the form 



Z = D9 t (t) exp 
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At T = 0, the (imaginary) time direction becomes infinite 
in extent, and the model is in the 3d xy universality class. 
The dualization procedure used for the classical 2d xy 
model can now be implemented on the 3d xy action as 
well. First we discretize the imaginary time direction 
in units of At and work on a three dimensional lattice. 
Introducing the variables m which only live on the spatial 
links, the action is 
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where c = a/ At. Rearranging terms as before we get 



where J — > J a 2 At, C — > Ca 2 / At and m — * m/a. The 
isotropic model is recovered when the discretization in 
time is such that At = \J C / J. Combining V x m 
(which is a scalar since m is two dimensional) and dm/dr 
into one three dimensional field we see that ean llOl is 
the partition function for the three dimensional loop gas 
model which is in the same universality class as the 3d 
xy models. A point to note here is that a vector field 
coupling via a three dimensional kernel does undergo a 
phase transition while a scalar field does not—. Alterna- 
tively, the loop gas model undergoes a 3d xy transition 
while the three dimensional Coloumb gas is always in a 
disordered phase. 

Alternate derivations^ map the model directly to the 
3d xy model instead of the loop gas model. We find 
that the procedure described above is preferable in the 
presence of dissipation. 
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In the dissipative 2d xy model, one adds to eqnJS] a 
term which generates in the equations of motion of the 
phase variable a term proportional to the velocity, 9. The 
dissipative part of the action is of the forro^ 
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where a — Rq/R where Rq = h/4e 2 . For a single 
Josephson junction the voltage drop is given by d{9\ — 
62) /dr. The power dissipated by a resistor connecting 
the two superconducting grains is {d{6\ — 9 2 ) / 'dr) 2 / R. 
Integrating over time one obtains total energy dissipated 
through the resistor. For an array of junctions the cor- 
responding term in the action is of the form in cqn llll 
Such a dissipative term can also be derived by assuming 
that the phase difference across every junction is coupled 
to an external bath of harmonic oscillators. Integrating 
out the bath the dissipative term is obtained in Fourier 
space as 

Sdiss =«£J]kri^ -0 3 »»| a . (12) 

The exponent 77 reflects the form of the spectral density 
I 



Z 



where we have redefined a — > aa 3 . The action in Eqnll3l 
has two possible phase transitions, depending on whether 
the capacitance C or the dissipation term oc a in the 
numerator of the second term in Eq. (|I3p dominates 
in long wavelength and low frequency limit. The former 
corresponds a critical point with the dynamic critical z = 
1, i.e. we recover the loop gas model which belongs to 
the 3d xy universality class. The latter corresponds, as 
we will show to z = 00, i.e. a fixed point with local 
criticality. The two terms have the same scaling form for 
z = 2. In the rest of the paper we analyze the dissipation 
driven transition 



The vector field m can be written as a sum of a lon- 
gitudinal, mi, and a transverse component, m t which by 
definition satisfy V x m; =0 and V -m t = 0. The vortic- 
ity field p v (r, r) is related to the transverse component 
alone and is given by 



I 

of the external oscillator bath at low frequencies. For a 
linear spectral density rj — 1 and a Ohmic dissipation 
is obtained^. It is indeed interesting to ask what other 
forms of dissipation is allowed. We will postpone the 
discussion of the effect of such dissipative terms to section 
X. 



II. EXCITATIONS OF THE QUANTUM 2D XY 
MODEL 

We can now carry through the program of introducing 
the vector field m for the quantum 2d xy model with 
Ohmic dissipation (77 = 1). The integration of the 9 
fields can be performed since the dissipative term is also 
quadratic. The effective model is now written as 
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V x m f (r, t) = p v (r, r) = p vc ^ 5 (r - 5(r- n) , 

(14) 

where (rj,Tj) are the location of the core of the vortices 
and p vc is the quantized charge of the vortex. 

Unlike the z = 1 theory where dm/dr is the two com- 
ponent current that acts as the additional topological 
defect, we find that the topological excitation in terms of 
which the problem with dissipation is most conveniently 
discussed involves only the longitudinal component and 
is given by a "warp" field p w (r, t) defined through a non- 
local relation in space, 
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J r — R 
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p w are events where the longitudinal component 
changes (see FigJTJ). The warp charge p wc is defined by 

Pw (R, r) = ^-p wc S(R - Ri)<J(T - n), (16) 

i 

in a way analogous to the definition of the vortex charge 
p vc above. We will show in Sec. (Ill) that the warp- 
charge is also quantized. 

Given Eq. (fT5|) . the relation between the Fourier- 
transformed warp field p w (k,ui n ) and the longitudinal 
component of the velocity field is 

• m ; (w„,t) = p w (k, u) n ) (17) 




FIG. 1: A "warp" is an event in time where the divergence of 
m field changes 

In terms of p v ,p w , the action, in the continuum limit, 
neatly splits into three parts: S = S v + S w + S' w 

S " = Z^£^M k "" )|2 ( 18 ) 

I 



x \p w (ku n )\ 

where 



O (K, U n ) - ^ + {c/c) ^ + a ^ fc2 (M) 

The factorization of the action in eqn !18l is only possi- 
ble due to the choice of variables made in eqnQTI With 
this choice, it follows from Eq. fTg)) there are logarith- 
mically interacting vortex charges in space with local 
interactions in time, and a logarithmically interacting 
warps in time with local interactions in space. The 
former is the standard Kosterlitz-Thouless-Berezinskii 
theor y 30 ! 31 while the latter is the charge representation 
of the Kondo problem 3 ^ near the Ferromagnet/ Anti- 
Ferromagnet quantum critical point. Quite unlike the 
z = 1 theory where Lorentz invariance is manifest, the 
space and time degrees of freedom are completely de- 
coupled in the theory describing the dissipation driven 
transition. Being orthogonal the charges p v and p w are 
uncoupled; the action is a product of the action over con- 
figurations of p v and of p w . Any physical correlations are 
determined by correlations of both charges. The third 
piece of the action S' w is nonsingular because it corre- 
sponds to the three dimensional Coloumb gas problem 
which is known to be disordered at all temperatures^ 3 -. 

The singularity induced by the presence of Ohmic dis- 
sipation is transparent in our formulation. Suppose, one 
works instead with the degrees of freedom appropriate to 
the 2+1 d xy model which are the vortex density p v and 
the vortex current, 3 V . The latter is defined as 



j v= 2x— (20) 
clt 

The action in eqnJT3] is 
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(21) 



The first two terms can be combined to give the dual 
representation of the 3dxy model. In the absence of dis- 
sipation, these terms lead to the superconductor to insu- 
lator phase transition. The third term introduces a new 
singularity in the vortex currents. Due to the continuity 
equation between p v and J„, see below, the two singu- 



larities cannot be decoupled as before and the analysis of 
the model is not straightforward. 

In terms of warps and vortices the vector field can be 
written as, 

m = p v z x ik/k 2 + p w ck/iu> (22) 
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and 



-p v uj\i/k 2 + p w z x ck 



(23) 



Separating the longitudinal and transverse components 
of the vector field m and the vector current J„ allows us 
to factorize the action. Notice that the new term in the 
vector current is transverse and the continuity equation. 
dp v /dr + V • J„ = is satisfied. 

Since warps are excitation that introduce longitudi- 
nal components in the vector field, alternate formulations 
could have been considered. In particular, we could have 
chosen to define monopoles (sources and sinks in two di- 
mensions) which are of the form 



HI = p rn 



(24) 



Such excitations determine the longitudinal component 
of the vector field and their relation to warps in momen- 
tum space is, 



p m (w,k) = — p w (w,k) . 



(25) 



In terms of monopoles one cannot separate the singular- 
ities and the simple description of the dissipation driven 
transition is no longer possible. The physics remains 
same as the scaling dimension of monopoles and warps is 
the same. We choose to work with warps as the analysis 
is more transparent and we do not have to consider the 
added complication of the renormalization of velocity. 



III. WARPS AS TOPOLOGICAL DEFECTS IN 
SPACE TIME 



corresponding 9 field has a discontinuity along a line go- 
ing from the core out to infinity. One can then proceed 
to the continuum to analyze the critical properties of the 
model . 

We see from the solution of Eq. (fT5|) that the field m/ 
due to a warp created at site R; and time Tj has the form: 



mj(r,r) = cp wc e(r - n 



Ri 



R,l 3 



(27) 



Analogous to the vortex, this solution is a quantized 
charge in a core which serves as the source of a radial 
velocity field falling off as 1/r 2 . To see this consider the 
solution to eq[T7]with m = V0, 



'(r,T) = e(r-Ti) ( 2nn 



a 



Ir-Ril 
(r-Ti) 0c, |r-Ri| < a 



v / , |r - Rj| > a 

(28) 



The boundary conditions chosen are 6(\r — Rj| = a) — 
6 C + 2irn and 6 (|r — Rj| = oo) = C . The first guaran- 
tees the quantization of the warp charge, imposed by the 
phase jump of 2im that occurs at |r — Rj| = a. In other 
words a warp is an event which produces a discrete phase 
jump across the core. The corresponding discontinuity in 
the phase is discrete, corresponding to a discrete field 
localized along a loop encircling the core. This localized 
discrete mi serves as the source of the radial field 
outside the core. As one moves along any radial direc- 
tion crossing this loop, which defines the boundary of the 
core, the phase jumps by 2nn. Thus the warp charge is 
discrete and satisfies the constraint that the vector field 
m is discrete. The second boundary condition where is 
chosen because we are only describing the changes with 
respect to the fields because of the warp event at r = . 



In this section we discuss the topological properties 
and the physical content of the warps. The definitions in 
Eqs. (fT7|) and (TT5|) say that the warps are source for the 
longitudinal component of the vector field m. Since m is 
an integer vector field, it is important also to show that 
this property is preserved by Eq. (fT5j) as well as that the 
warp-charge p wc is quantized. 

Let us first review the familiar case of the vortex in 
Eqnj26] In real space the vector field m t due to a vortex 
at R,- is 



m* (r,r) 



z x (r- R;) 
; Ir-Rd 2 



(26) 



Given the integer vector field m the vortex is defined by 
its quantized charge p vc and the singular azimuthal vec- 
tor field oc 1/r. This vector field is undefined inside a core 
radius a where the source of the velocity field with dis- 
crete charge p vc sits. The discreteness of m t is respected 
once such a singularity is identified. In particular, the 



DISSIPATION DRIVEN PHASE 
TRANSITION 



IV. 



The action for the warps, Eq. (Tl8j) . has been studied 
before in the context of quantum coherence of two state 
system coupled to dissipative bath a 37 i 38 . Let us introduce 
a core-energy A for the p w 's just as is done to control 
the fugacity of p v 's, the vortices. Next consider how the 
renormalization of a and A proceeds.—. Including the 
core-energy, the action S w is 



TZ~T7~\ \P wi ( w ™)| 2 + / dTA\p wi (T)\ 2 



Air w„ 



The RG equations for eqn[29]are well known 3 ^^^ 
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(30) 
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where A = Ar c and r c w (2J2-B) 1 / 2 is the short time 
cutoff. The critical point of interest is at a c = 1, where A 
scales to 0; for a < 1 the charges p w freely proliferate as 
"screening" due to A becomes effective, a > 1 represents 
the ordered or confined region. We are interested here 
only in the region a < 1. Well in the quantum critical 
region, the (singular part of the) propagator for p w is 



(p w (U3 n ) p w {-UJn)} = (1/47T \uj n \ T c y 



(31) 



Note that the form of the propagator in eqnJ3T] is valid 
only at the critical point. The crossover to the quantum- 
disordered or screened state is given for T = when 
u> = uj x which is of the order of the inverse of the charac- 
teristic screening time, which may be estimated similarly 
to Kosterlitz's estimate^ of the screening length in the 
xy-problem: 



■ 1 exp(-6/v / T — 



(32) 



where b is a numerical constant of 0(1). At finite tem- 
peratures and low frequencies, the crossover temperature 
T x is of 0(u> x ). On the ordered side a > a c , we have to 
consider the effect of the charging energy and the flow 
to the 2 = 1 critical point. In figd] we plot the phase 
diagram in the T — a plane where only the the region 
a > a c is correctly described. The crossover line is also 
shown where well above T x we expect to observe critical 
behavior controlled by the quantum critical point while 
below T x the system is phase disordered and hence an 
insulator. 



models. We are interested in the expectation value of 
exp(id ifl ) given by 



C. 



-idi 



Jd[9], 



C = 



fd[6]e- 



2m V 



m 



(33) 



(34) 



Notice that the exponential of the sum along the path is 
always one as m is an integer fields. We need to include it 
in our calculations to reproduce the correlation functions 
in the absence of dissipation. Since m is two dimensional 
vector, the contribution is only from the space like seg- 
ment of the path chosen. We can now proceed as before. 
Going over to imaginary time implies that the factor of 
1 is absorbed in the definition of the path integral. 



FIG. 3: Path between ip and jv used to compute the corre- 
lation function. 



To compute the correlation function we introduce a field 
if ^ that lives on the sites of the space time lattice. Its x 
and y component are 1 or depending on whether or not 
the path includes the link {xyp) to (x + ly/x) or (xy + lfx) 
respectively. Similarly we introduce a field ?7? which is 1 
or depending on whether or not the path includes the 
link (i/i) to (ip+ 1). The path we choose to compute 
the correlation function is shown in figO We can now 
integrate out the phase degrees of freedom as before to 
obtain the functional in terms of the vector fields alone. 



1/ot 



FIG. 2: Dissipative driven proliferation of warps occurs for 
1/a > 1. The power law scaling at criticality is cut off at an 
energy scale T x as discussed in the text. 



V. CORRELATION FUNCTION 

We can compute the correlation function for the order 
parameter field following a procedure which is a general- 
ization of the method developed for the classical 2DXY 



S — iS^j T 



+ 



L 2 /3 



(2tt) JG(k z , W „) [(-iuj n rf* /2) (ik t ■ m ln 



+ (-iki x if*) ■ (iki x m ; „) + c.c. 
+ {{C/c)loI + a\oj n \ k 2 ) (m • if) 

+ \u n ri /2 + k l -Jf*/2\ 



(35) 



Given the definition of vortex and warps, we can write 
the vector field as 



jkxz 



Pv 



k 2 



Pw k 

HO 



(36) 
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Replacing in the expression above we see that the vortex 
part of the correlation simplifies just as it did in the case 
of the action itself. 



space like component of the path) while the second term 
couples to time derivatives which are finite at (r',r) and 
(r',r'). 



S (Sal 



S w 

In 



-iuj n rp* /2) —cp u 



p w (k • If ) + c.c. 



({C) + cak 2 / \uj\) 

h 

■|w n »7 /2 + ki--^*/2| 
f 2ttJ . , 

S -p- (-»kj x 77 ) • p v + cc. 



The net correlation is a product of three distinct contri- 
butions. The term independent of vortices and warps is 
the spin wave contribution. Notice that the associated 
propagator is three dimensional and as such is non sin- 
gular. This reproduces the well known result that spin 
waves do not disorder the long range order in three di- 
mensions. The other two contributions come from vor- 
tices (the last line in eqn l37|) and warps. At the dissi- 
pative critical point, the singular fluctuations are asso- 
ciated with the warps. Thus the vortices and spin wave 
contributions are completely regular. To compute the 
contribution from the warps we follow a procedure sim- 
ilar to Jose et al. (— ). Consider the contribution from 
the warps, 



S w = S w + -±- y)(27r) JG(k,,^„) 



(-icj n <n°*/2) — Pu 



S w + 





ck 


1 dRdn 


— Pw 

U) 



(Rn) 



{G(rV,Rn)-G(ry,Rn)} 

(C + ack 2 /\uj n \) 



(37) 



dRdri 

h 

{G (rr, Rn) — G (r'r, Rn)} 
Fourier transforming back we get 



(Rn) 



(40) 



S — S v 



1 

Up 



(41) 



4 k r (r 3 )+jw Il (r„) 



iki-(rj-)+iw n (T M ) 



(C + acfc 2 /K|) 



( ki-(r 4 )+iw„(T M ) 



Since the action S w is quadratic, we can compute the 
correlation to give 



C(r-r',T- 

Pc 2 
F = - — 



G 2 (k,w n ) 



(C + ack 2 1 'H) 
ki 



p w (k • if ) + c.c. 



- t') oc exp(— F), 

dk (p w (w„) p w (-w„)) x 

n J 

k 2 

— (1 - cos (w„ (r - t'))) 



(42) 



(38) 



In real space the terms of the action can be rewritten 



'((C/c)+afe 2 /K|)' 



'((C/c) + afc 2 /K|)' 



(1 -cos(k/ ■ (rj-rj-))) 



ik r (ri 



+ e 



laj„(T M -T„) 



5,, 



dRdR'drdr' [ft-lf (R, r) 



(39) 



G ("Rr; R't'1 



(G + acfc 2 /|w|) 



fR'i 



ck 


(RV)" 


— Pw 


UJ 





9 T? 7 (R, r) G (Rr; R'r') 



where the derivatives are taken along the path and the 
term [..] (R, r) is meant to represent the fourier transform 
of the terms within the square brackets. Consider the 
derivatives along the path of the fields If and 77 . They 
represent the change in value of the fields computed from 
one point to the next along the chosen path. Since the 
first term in eqn 1391 couples only to the spatial derivative 
it has finite values only at (r, r) and (r',r) (ends of the 



ik,-(i\-r i )+tu< n (T„-T 1 ,) _ ! 



The order parameter correlation function has three con- 
tributions from the warps. The first term in the square 
brackets in eqn[42] is from the time like segment of the 
path; the second term corresponds to the space like seg- 
ment of the path and the last term is an interference term 
between the two. At the critical point the correlation of 
the warps is given by eqn l31l Consider the space like seg- 
ment. It is singular (of the form (1/ |w n |)) as oj goes to 
zero. This implies F — > 00 and the correlation function 
is zero unless = rj. Thus the correlation function is 
local. For = rj the last two terms are zero and only 
the time like path contributes to the correlation function. 
Performing the integral over momentum k, the correla- 
tion function is given by G(0, t—t') oc exp(— F(0, t—t')), 
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where 

F(0, r-r<) = -2nT £ ^sK^)) ^ ( ^ 

(43) 

To perform the matsubara sum we consider the analytic 
properties of the logarithm in eqn l43l The branch cut 
introduced on the real axis implies that on integrating 
over a contour infinitesimally above the real axis from 
— oo to oo and infinitesimally below the real axis from 
oo to — oo the real part cancel while the imaginary part 
enforces a sum over positive frequencies giving 

Such correlation functions have been calculated in other 
contexts. In particular using the calculation by Ghaemi 
et al42, the spectral function of the correlation is 

Zmx(q,w) = f c tanh(w/2r), u < t" 1 . (45) 

f c is 0(t c ). The order parameter correlation at a = 1 
is local in space and power law in time. Thus the dis- 
sipative 2DXY model describes a new type of quantum 
criticality where the correlations are local in space and 
are functions of ui/T, with a cut-off provided above. Pre- 
cisely such criticality had been postulated to describe the 
"Strange metal" or marginal fermi- liquid region of the 
phase diagram of high temperature superconductors^. 

VI. ZERO TEMPERATURE PHASE DIAGRAM 

So far we have concentrated on the dissipation driven 
transition. We now consider the phase diagram in the 
JC — a plane. Our analysis clearly shows that the physics 
at finite dissipation and at zero dissipation cannot be con- 
tinuously connected. This is in accord with the sugges- 
tion in Ref.(lg), that the zero dissipation limit, described 
by the 3d—xy model at zero temperature, is in a different 
universality class compared to the model at infinitesimal 
dissipation. The other limit where the physics is easy to 
analyze is at infinite dissipation. In this case, the fluctu- 
ations in time are energetically unfavored and the model 
maps to a 2D — xy model. We now make a guess on the 
structure of the rest of the phase diagram based on the 
well-defined theory in the limits discussed. 

In reference [6] it is argued that a weak coupling analy- 
sis shows that the entire Josephson coupling renormalizes 
to zero at small dissipation. Furthermore a self consis- 
tent calculation allows one to map the phase boundary 
within this scheme. Such an analysis does not capture 
the physics of topological defects that are necessary to 
determine the entire phase diagram. Given our identifi- 
cation of the relevant degrees of freedom and the fugac- 
ities associated with the new excitations, we propose an 



alternative understanding of the physics of the dissipative 
2d — xy model. 

We have explicitly shown above that the critical theory 
has two independent singularities. Thus the true phase 
diagram can be determined by studying the flow of stiff- 
ness, fugacity of vortices and fugacity of warps. The only 
algebraically ordered phase is one in which both the vor- 
tices and warps are bound. This occurs for 1/a < 1 and 
JC > 2/ir. Furthermore this phase is continuously con- 
nected to the 2d XY phase at infinite dissipation. Thus 
one expects a phase boundary to exist that starts at 
JC = oo and a — 1 and approaches JC = and a = oo 
asymptotically as shown in figlU For 1/a > 1 the fugac- 
ity of warps grows and the stiffness vanishes. Similarly, 
for JC < 2/ir the vortices proliferate driving the stiffness 
to zero. We find no evidence for the ordered phase for 
small JC and a > 1 as in reference [6] where the large 
fugacity of the vortices was ignored. Similarly, for large 
JC a minimum value of a is required to overcome the 
fugacity of the warps. In other words for any given value 
of JC as one increases a the fugacity of warps, which 
is infinite at infinitesimal a, remains large until we cross 
the phase boundary, upon which it renormalizes to zero. 
Whether the resulting phase is algebraically ordered or 
not depends on the fugacity of vortices. 



♦ 




1/a 

FIG. 4: The zero temperature phase diagram in the 1/a — JC 
plane for the dissipative 2DXY model. For a = 0, the sys- 
tem undergoes a superconductor to insulator phase transition 
with a 3DXY critical point. For infinite a, the dynamics is 
quenched and the we recover the 2DXY behavior. There are 
two phases for finite dissipation. The disordered phase is char- 
acterized with a diverging fugacity of warps and vortices. The 
algebraically ordered phase is one in which both the vortices 
and warps are bound, and the correlations are power laws 
separately in space and time with continuously varying expo- 
nents. 



Consider the regime where warps proliferate. The van- 
ishing stiffness implies that such excitations should also 
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lead to screening in the interaction of the vortices. In 
analogy with the disordered phase in the KT theory, we 
should introduce a mass for both vortices and warps. 
This in turn will lead to unbinding of vortices as well. 
We believe that for small JC, this proliferation of warps 
and resultant screening of interaction among vortices will 
push the phase boundary away from a = 1 to larger val- 
ues of a. At large JC we expect our theory to be valid 
and a dissipation driven transition will occur at a = 1. 
A schematic phase diagram in shown in figHl 

At first sight it might appear that a long range interac- 
tion in time of the form 1/ |r — t'| 2 is inconsistent with 
a disordered ground state, where the correlation func- 
tion decays faster than the interaction. The reason such 
a state is allowed is the fact that the existence of ex- 
citations which correspond to phase slip events in time 
introduces a second time scale in the problem. The com- 
petition between the short time scale physics associated 
with the fugacity of the warps and the long range interac- 
tion between them is precisely the reason for the observed 
phase transition. 



VII. ASHKIN TELLER MODEL AND THE 
CUPRATES 

In reference [1] we have argued that the long wave- 
length theory of a model for the Pseudogap state for the 
Cuprate superconductors is the Ashkin- Teller model, The 
quantum version of this model with dissipation is in the 
universality class of dissipative quantum 2d — xy model. 
The Marginal Fermi Liquid 24 properties of the optimally 
doped system originate from the fluctuation spectrum 
of the local quantum critical point. In this model the 
phase variable 9 represents the possible current loop or- 
der within a given unit cell and the four fold anisotropy 
restricts it to four values representing states that break 
time reversal and one of two possible reflection symme- 
tries. The origin of the dissipation in this system is due to 
the coupling of the order parameter, which corresponds 
to the coherent part of the current operator, to the in- 
coherent part of the current operator which correspond 
to fermions near the Fermi surface. For long-wavelength 
fluctuations such a coupling leads to Ohmic dissipation of 
the form in eqn llli while term of the form in eqn !56f see 
below) is not generated. Detailed derivation of the long 
wavelength theory will be provided in a later publication. 

The classical Ashkin Teller model is defined in terms of 
two Ising degrees of freedom per unit cell. If we define oi 
and Ti as Ising variables taking ±1 values to parameterize 
the currents along the the x and y direction, the four 
resulting states in each unit-cell (see figlS]) form the basis 
of the low energy sector of the long wavelength theory. 
The Hamiltonian of the classical Ashkin Teller model is 22 

Hat = - ^ t^ 2 + Ti7 j) + J^crjTiTj] (46) 

<ij> 



+ 1 .1 



a 




o ► • -o o ►o 

I. \i 




• Cu o O 

FIG. 5: The four domains of the circulating current phase 
is shown. They can be parameterized by two Ising variables 
which represent the currents in the x and y directions. 



For J4 = the model reduces to two decoupled Ising 
models which is known to undergo a Gaussian phase 
transition 4 ^. The analysis has been extended to map 
out the entire phase diagram by observing that J4 is a 
marginal operator in this theory 4 ^. For —1 < J4/J2 < 1 
the system undergoes a Gaussian phase transition form a 
ferromagnetically ordered state to a paramagnetic state 
with continuously varying exponents. Quite remarkable, 
for —1 < J4/J2 < 0, the specific heat exponent a is 
negative. The Ashkin Teller model explains why no 
singularity is observed in specific heat measurements at 
the pseudogap temperature accompanying the breaking 
of time reversal symmetry. Defining cr,; = cos (0i) and 
Ti = sin (9i) we can rewrite the Hamiltonian as 

H AT = - t 2j2 cos ~ e i) + Ji cos _ 29 i^ 

<ij> 

(47) 

where OiS take values 7r/4, 37r/4, 5tt/A and 77r/4. The 
constraint can be implemented as a four fold anisotropy 
term, 

H AT = - t 2 1/2 cos + J ± cos ( 2(? * _ 29 j)] 

<ij> 

+ ^/i 4 cos (40*) (48) 

i 

Given this mapping of the Ashkin Teller model to a xy 
model with four fold anisotropy, we now analyze the effect 
of hi on the local quantum critical point of the dissipative 
2dxy model. 
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VIII. EFFECT OF FOUR FOLD ANISOTROPY 

The discrete nature of the underlying degrees of free- 
dom of the Ashkin Teller model enforce an Ising type or- 
der at low temperatures. On increasing temperature, it 
undergoes a Gaussian transition into a state which can be 
described by proliferation of vortices. It is a unique prop- 
erty of four fold anisotropy that the disordering, which 
occurs via proliferation of vortices, and ordering, which 
is enforced by a diverging anisotropy field, happen at the 
same temperature 2 ^. We now analyze the effect of such 
an anisotropy field on the local quantum critical point. 
To do so we introduce the term H4 = J^i hi cos (46^) in 
the Hamiltonian. To handle such a term in the action we 
employ the following approximation scheme 

Pi 

where pi is an integer field that lives at each site and 
1J4 = hi/2. For large values of j/4 the approximation is 
reliable as the sum will be dominated by the terms with 
Pi = 0, ±1. Eq|49lintroduces two new terms in the action, 
one linear in the phase and the other independent of it. 
We can proceed as before in integrating out the phase 
degrees of freedom to get a new action which is of the 
form 

1 x / k 

S = S V + S W + — )Gc 2 47r 2 J-p U) (k,w)p(-k,-w) 

+ b(k,^)| 2 ) +ln(y 4 )|p(k, W )| 2 . (50) 

Given that the coupling between the integer field p and 
the phase variable is linear we can follow a self consis- 
tent procedure to calculate its effect on the dissipative 
fixed point. In this context it is useful to recall that in 
the absence of dissipation, the disordered side is char- 
acterized by the anisotropy field going to zero, i.e. be- 
ing marginally irrelevant. Furthermore, in the absence 
of anisotropy the dissipative transition leads to a state 
which is characterized by the proliferation of warps. In 
other words, the most singular part of the action is the 
interaction between warps in time. Since the action for 
the integer field, p, is quadratic, the cumulant expan- 
sion implies that the a new term in the action for warps 
appears of the form 

~ 1 / k \ ^ 

54 = 2L^3^ ( 47r2jGc2 ^ ) P^{K^)Pw{-K-^)^< 
(p(k,w)p(-k,-w)>. (51) 

The correlation is obtained from a self consistent solution 
of the coupled equations shown diagrammatically repre- 
sented in Fig (HI 

A quadratic equation is obtained for each of the prop- 
agators. Since we are looking at the effect of the 



= + 

(a) 

= + 

(b) 

FIG. 6: Equations for he full Greens function. The thick 
lines represent the full propagator {p (k, uj) p (— k, — uj)) and 
the thick dashed lines represent (p w (k, uj) p w (— k, — uj)), while 
the thin full and dashed lines represent the same propagators 
computed in the absence of the linear coupling with the warps. 

anisotropy on the effective action for the warps, we can 
solve these equations to obtain the most singular con- 
tribution to the propagator in Eqnf5U Defining G p = 
(p(k,w)p(-k, -ui)) and G w = (p w (k,u)p w (-k, -uj)), 
symbolically the equation satisfied by G p can be written 
as 

A 2 GjJ,G 2 — Gp + Gp = 0. (52) 

where G° and G° are the corresponding propagators 
evaluated in the absence of the linear coupling and 
A = (47r 2 JGc 2 k/oj). At the quantum critical point, the 
most singular contribution for the warps is logarithmic 
in frequency, while the action for the integer fields is not 
singular. This is again due to the fact that the integer 
fields interact via Colomb Kernel in three dimensions. 
Thus the most singular part of the solution is of the form 
G p w 1/A 2 [G°]~\ Substituting in EqnElJ the effect of 
the anisotropy is to introduce a term in the action of the 
form 

^ w 2^g £ t G ° ] ~' Pw ( k ' ^ Pw (~ k ' ~ u ) ■ ^ 

Thus all anisotropy does is to modify the coupling con- 
stants but the local quantum criticality remains and the 
form of the correlations at the critical point are still local 
in space and power law in time. 

IX. ORIGIN OF DISSIPATION 

In addition to the bosonic modes associated with the 
order parameter, the low energy sector of the cuprates 
has fermionic degrees of freedom associated with elec- 
tron near the fermi surface. Quantum fluctuations of the 
order parameter correspond to flipping between states of 
the Ashkin- Teller model. From the mapping above of 
the Ashkin Teller model to an xy model with four fold 
anisotropy, we note that a flip operator is the quenched 
angular momentum operator. For the cuprates, the anal- 
ogous operator can be shown to correspond to the curl 
of the fermionic current^. The coupling between the 
fermionic degrees of freedom and the bosonic critical 
modes is of the form 
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H coup ~ U(r) (V x j) (54) 

where U is an operator which flips the state of the Ashkin 
Teller model (analogous to d/d6 in the continuous limit) 
and j is the fermionic current. Integrating out fermions 
leads to an Ohmic dissipative term in the action for the 
order parameter. The reason is that the current current 
correlation is the conductivity which goes as \uj\. Un- 
like Josephson junction arrays, the origin of dissipation 
in cuprates is the presence of both fermionic and bosonic 
degrees of freedom in the low energy long- wavelength the- 
ory 



X. COUPLING TO OTHER DISSIPATIVE 
BATHS 

The form of dissipation considered so far is physically 
motivated to reproduce a term linear in the gradient of 
the phase variable in its equation of motion. Furthermore 
it can also be understood as arising from dissipation via 
resistors connecting the islands on which the phase de- 
grees of freedom live. In principle one can also study the 
effect of other forms of dissipation that are allowed by 
symmetry. In particular, we can introduce a term which 
leads to suppression of density in each grain instead of 
the phase difference between nearest neighbor grains as 
in eqn llll 




To understand the effect of such a term on the dis- 
sipation driven critical point we consider the effect of 
integrating out frequency shells as before. We find that 
the coupling constant j n scales to zero since it has a di- 
mension of —2. In other words such a form of dissipation 
is irrelevant near the local quantum critical point, while 
for a z — 1 critical point such a term is marginal. 

A dissipative mechanism where the external bath di- 
rectly couples to the phase variable leads to an action of 
the form 

i J \t-t'\ 

The coupling is marginal at the local quantum critical 
point. A study of the critical properties showed the ex- 
istence of a line of fixed points controlling the transition 
from the superconducting to the insulating stated, but 
the approximation made were called into question^. A 
weak coupling study of the model has been reported 51 
and a strong coupling analysis will be the subject of fu- 
ture investigation 



XI. RELATIONSHIP TO OTHER MODELS 

With the identification of the correct degrees of free- 
dom we have been able to access the dissipation driven 
critical point that has been previously shown to exist 
within a perturbative analysis. Furthermore, fourfold 
anisotropy does not alter the local criticality obtained 
at this quantum critical point. Thus the dissipative xy 
model defines a new universality class. Classical models 
which undergo a Gaussian phase transition from an or- 
dered to a disordered phase will belong to this universal- 
ity class provided their quantum generalization include 
Ohmic dissipation. In this context we can look at the 
physics of six and eight vertex models. 

The phase diagram of the six vertex model has two 
ordered phases, the orbital ferromagnet and the orbital 
antifcrromagnet, and a power law phase. The high tem- 
perature phase is the power law phase and at low temper- 
atures undergoes a transition to one of the two ordered 
phases. It is known that there exists an essential singu- 
larity at the phase transition and the physics is similar to 
the xy model except that is inverted. The high temper- 
ature phase is power law correlated while the low tem- 
perature phase has finite correlations^. Unlike the dissi- 
pative XY model considered here, the inverted character 
of the phase transition, while allowing for a phase transi- 
tion without a singularity in specific heat*£, is unlikely to 
belong to the universality class of local criticality. How- 
ever, the low temperature physics of this model including 
quantum generalization allowing for sources and sink is 
yet to be analyzed. The classical eight vertex model on 
the other hand has in its parameter space a high temper- 
ature disordered phase which can undergo a KT transi- 
tion to an ordered phased. Such models also share the 
feature that the specific heat is nonsingular at the phase 
transition. Quantum generalization of this model with 
dissipation have not been studied. In particular it is not 
clear which operators in the classical theory are equiva- 
lent to the warps defined in the continuum. A detailed 
analysis of the conditions for classical eight vertex mod- 
els to exhibit local criticality is beyond the scope of this 
paper. 

Local quantum criticality closely related to the form 
proposed for the Cuprates^ and derived here has also 
been used to fit the spectral function near the quantum- 
critical point of a heavy-fermion^ measured by neutron 
scattering. Some theoretical calculations^ emphasizing 
the deconstruction of the single impurity Kondo effect 
near the quantum critical point^ have been performed. 
A general point based on the present work appears worth 
making: Local quantum criticality requires showing that 
near the quantum-critical point, there must be variables 
in terms of which the action for low energies is separable 
exactly into orthogonal parts which interact only in space 
without retardation in time, and those which are power 
law in time but interact only locally in space. It ap- 
pears to us that only dissipation driven quantum critical 
points can have this property. At this point, it does not 
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appear clear to us in terms of which variables, the gen- 
eral Action of the heavy fermion problem with dissipa- 
tion, local Kondo interactions, interaction between local 
moments, magnetic order with attendant deconstruction 
of the Kondo effect, etc., may be cast into the simple 
structure due to which we have been able to solve the 
problems addressed in this paper. 

XII. SUMMARY 

In this paper we have provided a detailed analysis of 
the dissipation driven phase transition in the quantum 2d 
xy model. There exists a critical point in the phase dia- 
gram where the correlation are shown to be local in space 



and power law in time. This is a new paradigm in criti- 
cal phenomena which does not appear in the universality 
classes of classical dynamical critical phenomenal. The 
quantum disordering due to the proliferation of a new 
of class of topological defects, the warps, which inter- 
act logarithmically in time but local in space allows such 
criticality. We have discussed that the results here are 
only asymptotically true and only valid for certain forms 
of dissipation. A verification of the results by quantum 
monte-carlo calculations is highly desirable. 

A principal outcome of the paper is to derive the phe- 
nomenological assumptions of the marginal fermi-liquid 
theory with which the critical properties of the Cuprates 
in the anomalous metallic regime have been understood. 
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